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ABSTRACT 

The idea that GRBs originate from uniform jets has been used to explain numerous observations of breaks in 
the GRB afterglow lightcurves. We explore the possibility that GRBs instead originate from a structured jet that 
may be quasi-universal, where the variation in the observed properties of GRBs is due to the variation in the 
observer viewing angle. We test how various models reproduce the jet data of Bloom, Frail, & Kulkarni (2003), 
which show a negative correlation between the isotropic energy output and the inferred jet opening angle (in a 
uniform jet configuration). We find, consistent with previous studies, that a power-law structure for the jet energy 
as a function of angle gives a good description. However, a Gaussian jet structure can also reproduce the data 
well, particularly if the parameters of the Gaussian are allowed some scatter. We place limits on the scatter of 
the parameters in both the Gaussian and power-law models needed to reproduce the data, and discuss how future 
observations will better distinguish between these models for the GRB jet structure. In particular, the Gaussian 
model predicts a turnover at small opening angles and in some cases a sharp cutoff at large angles, the former 
of which may already have been observed. We also discuss the predictions each model makes for the observed 
luminosity function of GRBs and compare these predictions with the existing data. 



1. INTRODUCTION 

One of the outstanding problems in the field of GRBs is un- 
derstanding the extent to which these events are beamed, as 
well as the structure or configuration of the jet that produces 
the burst and its afterglow. The simplest model is a uniform jet, 
in which it is assumed that all parameters (density n, Lorentz 
factor r, magnetic and electron equipartition factors e# and 
e e , etc.) are constant throughout the jet. This jet can be de- 
scribed by an opening angle Oj (or alternatively a solid angle 
2ir(l - cosflj) ~ irO 7 -,). Under this assumption of uniformity 
throughout the jet, Frail et al. (2001) inferred the jet angle to 
a sample of GRBs with observed afterglows, from a break ob- 
served in the afterglow light curve. In a uniform jet model, this 
break occurs about the time when the GRB ejecta has slowed 
down enough so that the relativistic beaming angle of the radi- 
ation, ~ 1/r, becomes greater than 6j (e.g. see Rhoads, 1997, 
Frail et al., 2001). Using the inferred jet angle, the measured 
flux and redshift to each burst, Frail et al. (2001), and then 
Bloom, Frail, & Kulkarni (2003; hereafter BFK) with a larger 
sample, were able to determine the emitted energy E of each 
burst. Remarkably, they found that the GRBs in their sample 
exhibited very little dispersion in E (see Figure 1 of BFK); in 
other words, they found that the isotropic equivalent energy of 
the GRB, Ej SO , and the jet opening angle, 9j, adhere to the re- 
lationship Ej S0 6 2 j ~ constant. This intriguing result lends some 
credence to the possibility that a uniform description of a GRB 
jet may be valid. 

However, an alternative model was suggested by Rossi et al. 
(2002) and Zhang & Meszaros (2002; hereafter ZM). They sug- 
gested that in fact GRBs may have a structured jet configura- 
tion, in which parameters such as the emitted energy can vary 
as a function of angle from the jet axis. All GRBs may then 
have approximately the same ("quasi-universal") jet profile but 
appear different because of different orientations of the GRB jet 
to the observer (or in other words, varying observer angle #,,). 
In this model, a break in the afterglow light curve is still ob- 



served at the time the Lorentz factor slows to a value T ~ 1/8 V 
(see Rossi et al., 2002, and ZM for more discussion). There are 
several advantages of this model. It has the appeal that it allows 
for uniformity among the GRB population, and also makes def- 
inite predictions about the distribution of observed break times 
in the GRB afterglow light curve (Perna, et al. 2003), as well as 
the observed GRB luminosity function (ZM). The uniform jet 
model has to require different bursts having different opening 
angles, and it has no predictive power for the observed distri- 
bution of these opening angles or the GRB luminosity function. 
In addition, realistic simulations (e.g. W. Zhang, Woosley & 
MacFadyen, 2002) naturally predict jet structure from the col- 
lapsar scenario, so it is essential to investigate the possible jet 
structures in the quasi-universal picture. Motivated by the Frail 
et al. result {Ej S0 8j =constant), the most straightforward univer- 
sal jet model would be a power-law with an index of k = -2, as 
Rossi et al. and ZM have discussed. However, ZM observed 
that in terms of interpreting the lightcurves and jet breaks, one 
does not need to abide by any power-law jet configuration. In 
particular, they suggested that a Gaussian structure (and more 
general configurations) may describe the GRB jet. It has also 
been shown that a simple power-law model violates some of 
the observational data. For example, Granot and Kumar (2003), 
found that this type of model cannot reproduce the observed af- 
terglow light curve, while Kumar and Granot (2003) have found 
that a Gaussian structure for the universal jet model does a bet- 
ter job. 

In this paper, we test how well the quasi-universal jet model 
can reproduce the existing data. In particular, we test power-law 
and Gaussian models for the emitted energy, e(9), as a function 
of angle from the jet axis. Because it is unphysical to expect that 
all GRBs have exactly the same jet structure, we allow for real- 
istic scatter in the parameters of the models and place limits on 
this necessary scatter. The paper is organized as follows: In §2, 
we describe the data and any possible selection effects that may 
play a role in the results. In §3, we show how pure power-law 
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and Gaussian models fit the data of BFK. In §4, we introduce 
scatter into the parameters of each model and show how much 
dispersion is needed in these parameters to accurately repro- 
duce the data. We also discuss predictions the Gaussian model 
makes in the Ei SO - Oj plane, which may be tested with future 
observations (and may have already been observed). In §5, we 
discuss what each model predicts for the observed luminosity 
function and how this compares to existing data. A summary 
and conclusions are presented in §6. 

2. DATA 

Our data is taken from BFK, which provides an isotropic en- 
ergy estimate for 27 bursts. Of these bursts, 16 have a clear 
break in the afterglow light curve from which a jet opening 
angle (in the uniform jet model) can be inferred. These data 
are shown in Figure 1 (square points). Of the 11 bursts elimi- 
nated from the original sample, 8 have upper (or lower) limits to 
9j, based on a steeper (shallower) light curve observed at later 
(earlier) times, but no break directly observed. The other three 
bursts were well described by a single power-law throughout 
the extent of the observations. As discussed in the introduc- 
tion, they found in this sample a very small dispersion for the 
energy emitted from a uniform jet of opening angle Oj. This 
translates to a negative correlation between the isotropic equiv- 
alent energy Ei SO and the jet opening angle, such that Ei so 2 ~ 
constant. We performed all of our analysis below on both the 
nominal 16 bursts as well as the sample where the 8 limit bursts 
are included (where we took the limit as the nominal value), 
which gives us some idea of how eliminating these bursts bi- 
ases the sample and our results. We find in fact the results do 
not change qualitatively or quantitatively in either case, so we 
consider our sample with the limits removed reasonably com- 
plete. We present results for the sample of 16 bursts for which 
a definite jet opening angle has been estimated. 

Of course, there could be additional underlying selection af- 
fects biasing the entire sample that must be considered. This is 
a concern, for example, if there is a selection against detecting 
high Eiso, high Oj bursts or low Ei SO , low Oj bursts. Such selec- 
tion effects could produce an artifical correlation in the E iso -0j 
plane. This latter bias could be particularly worrisome. First, 
a low Oj implies an early jet break time, and a steepening of 
the light curve that may lead to a missed afterglow detection. 
More importantly, however, is that a low E iso may be missed 
due to the flux limit of the detector. For a burst at some red- 
shift z, the limit on E iso is given by E iso j im = Fi im 4-n d 2 (1 + z), 
where Fu m is the limiting fluence (time integrated flux) of the 
detector and d is the metric distance as a function of red- 
shift. For bursts at a redshift of 1, for example, we find that 
Eisojim ~ 10 50 (F; (m /3 x 10" 8 ergcm" 2 ) erg. This is well below 
the values of E iso in our sample. [We also point out that there 
is some evidence that bursts have higher Ei SO at higher redshifts 
(Lloyd-Ronning, et al. 2002, Amati et al., 2002, Yonetoku et al. 
2003, Graziani, et al. 2003), which helps reduce the severity of 
the flux selection.] In addition, the expression for computing Oj 
(see equation 1 of Frail et al., 2001) is very weakly dependent 
(to the 1/8 power) on Ej S0 . This means that the functional form 
of the truncation in the E iso -6j plane goes as E iso oc 6*J 8 , which 
is far steeper than the correlation observed. As a result, we 
conclude that selection effects are probably not producing the 
correlation between E iso and Oj in the BFK data, and that this 
is a reflection of a real physical effect in GRBs. Under the as- 
sumption that this correlation is real, we explore configurations 
of the energy as a function of angle e(0) in the quasi-universal 



configuration. 

3. QUASI-UNIVERSAL JET MODELS 

In the universaljet model, the oberver angle V takes the place 
of the jet opening angle Oj. Then - as described in ZM - Ei SO (0j) 
in the uniform jet model translates to e(0), the energy as a func- 
tion of angle from the jet axis, in the quasi-universal picture. In 
principle, we can then fit Ei SO (0j) with various models to con- 
strain the possible quasi-universal jet structure e{0). As dis- 
cussed in the introduction and in ZM, two plausible models for 
the jet structure e(0) are a power-law and a Gaussian. We de- 
scribe each in turn below. 

3.1. Power-law 

One possible model for the structure of a universal jet is a 
power-law: 

e(0) = e o {0/0*y\ (1) 

for > C , where C is some cutoff to prevent divergence (see 
also ZM). The solid line in Figure 1 shows the best power-law 
fit to E Lw (0j) oc e(0). The best fit gives an index k of -1.9 ± .1. 
We emphasize that the reduced x 2 of this fit is unacceptably 
high (<~ 8) because of the intrinsic scatter in the data (we will 
return to this point when we compare Gaussian models to the 
data). Nonetheless, a simple power-law structure describes the 
general trend of the data well. Of course this result is not sur- 
prising since the Frail et al. (2001) and BFK result shows that 
EisoO 2 ~ constant (or E iso oc 0~ 2 ). Although this model appears 
to describe the data, it is worth investigating what other types of 
jet structure may produce the BFK data. In particular, Granot & 
Kumar (2003) ran simulations which showed that a power-law 
jet structure cannot reproduce the afterglow light curves - there 
is an unobserved flattening both right before the jet break and 
also at late times (after the jet break). A Gaussian, model, on 
the other hand, does reproduce the features observed in the af- 
terglow lightcurves (Kumar & Granot, 2003). We explore this 
model in the next section. 

3.2. Gaussian 

A Gaussian model, 

e(0) = e o e- s2 ^> (2) 

is also a physically reasonable suggestion for the structure of a 
universaljet. In this model, the energy peaks at = (i.e. down 
the center of the jet axis) and o can be considered a character- 
istic width of the jet. The dotted line in Figure 1 shows the best 
fit Gaussian jet structure, with a characteristic jet width o of 
13 degrees. By eye, one can see that this is an unacceptable de- 
scription of the data and indeed this is confirmed formally (with 
reduced \ 2 ~ 20). [We note that when a Gaussian + power-law 
model is tried, the Gaussian component is made negligible in 
the fit routine, consistent with the above results.] 

However, we point out that it is unphysical to expect that a 
jet's energy structure (as a function of angle from the axis) will 
be exactly the same from GRB to GRB. A much more natural 
and physical model is to allow a quasi-universal configuration 
- i.e. to allow some scatter in the parameters of each of these 
models. It is possible, then, that more general configurations 
(other than a simple power-law) can better describe the data. 
We explore this possibility in the following section. 



3 



3.3. Varying the Model Parameters 

We investigate how the Gaussian and power-law models im- 
prove in describing the data when we allow the parameters 
(such as the power-law index or mean of the Gaussian distribu- 
tion) to vary. We put limits on how much the parameters in each 
model much vary to reproduce the scatter in the observed data. 
As we will show below and as is consistent with our results in 
the previous section, a power-law model does a good job of de- 
scribing the data and requires only small intrinsic scatter in the 
model parameters. However, the Gaussian model can also re- 
produce the observed data well; futhermore, in most cases this 
model predicts a sharp cutoff in the Ei so - 9j plane at large 9j as 
well as a turnover at very low values of 9j, the latter of which 
may already have been observed. We discuss this in more detail 
below. 

3.3.1. Method 

To test this scenario, we allow one parameter (or two when 
noted) of each model to vary according to either a log-normal 
or normal distribution. For example we can allow the normal- 
ization of the Gaussian model e a to vary as a log-normal dis- 
tribution rather than be a single fixed value. For each value of 
9, 1 we draw the normalization e a (in this example) from a log- 
normal distribution with some mean and scatter, and keep the 
mean of the jet width 9„ fixed to some value. We then compute 
the energy using this normalization and o according to equa- 
tion 2. Because the normalization is drawn from a distribution 
with some scatter, this will produce some scatter in the energy 
as a function of angle. The point is to find out how much scat- 
ter is needed in each parameter of the models to reproduce the 
data. Besides requiring the data to qualitatively look similar, 
by reproducing the energy range and scatter of the BFK data, 
for the Gaussian jet model we would like the simulated data 
to be fit by a power-law with an index consistent with the fit 
to the BFK data (~ -2), with a \ 2 ~ 4- 10 (similar to that of 
the data). Finally, we run a Kolmogorov-Smirnov (KS) test to 
check whether the two distributions (data vs. simulations) can 
be considered statistically similar. 

We comment that for all parameters but the index k in the 
power-law model, the log-normal distribution is the probably 
the most physically reasonable distribution for adding disper- 
sion to the parameters of each model. It accurately reflects 
the dynamic range of variation seen in the data and avoids un- 
physical negative values for the energy that can arise when a 
simple normal distribution is used for the parameter variations. 
Also, if GRBs from progenitors over a wide mass range, e.g. 
10-300Mq, a log-normal distribution is more natural to de- 
scribe the dispersion in the physical parameters. If, however, 
they arise from progenitors of a very narrow mass range, a 
normal distribution for parameter variations may be appropri- 
ate. Although (except for the power-law index k) we consider 
the log-normal variations to be physically more reasonable, we 
present the results for both distributions below; our figures, 
however, show only the results of the log-normal variations (ex- 
cept, again, for the power-law index k). 

3.3.2. Results 

1. Power law: 

First, we explore how much scatter is necessary in the power- 
law model to reproduce the scatter of the BFK data. We 

'Our 9 values are taken from the BFK data directly. However, we find quali 
distribution in linear space. 



consider the case where the power-law index k in equation 
1 is not a constant but derived from a normal distribution 
oc zxp~ (k ~ <k>) l la i< . We choose a mean < k > of -2 and explore 
various cr* values that can reproduce the data. For a 0* = 1° 
in equation 1, a scatter of o> ~ 0.7 adequately reproduces the 
BFK data. This is shown by the circles in Figure 2 (the BFK 
data are the square points). The bottom line is that the scatter 
in the power-law index need be on the order of ~ 0.7 to re- 
produce the scatter in the BFK data. We also vary the normal- 
ization e a as a log-normal cx e - 1 exp" (lg(e " ) " <lg(e " )>) / 2<T w«o). We 
find that a mean of < lg(e o /10 50 erg) >= 4.8 and a dispersion 
oigfe/io^erg) = 0-6 adequately reproduces the data (in the sense 
of reproducing the scatter and preserving the k = -2 power-law 
behavior) as shown by the circles in Figure 3. When we vary 

e a according to a normal distribution oc exp - ^" -- * 6 "^ / 2<7e », we 
find that < e o /10 50 erg >= 1.6 x 10 4 and a €a =< e > does an 
adequate job of reproducing the scatter. In all cases, we have 
taken 0* in equation 1 equal to 1°. 
2. Gaussian: 

In this case, we tried varying each of the parameters, e„ and 9 
of the Gaussian model in equation 2, according to a log-normal 
and a normal distribution. For example, as described below, we 
let the dispersion or characteristic width of the jet 9 be derived 
from a log-normal function with a mean < lg(6>„) > and stan- 
dard deviation cri g (e„), while e is fixed to some constant value. 
In addition, we try a model in which the normalization e is 
taken from a log-normal distribution and then 8„ varies in a 
correlated way according to the equation e 9 2 ~ constant. This 
is motivated by the fact that the Frail et al. (2001) and BFK 
data imply that the total GRB energy is approximately constant 
(in a uniform jet configuration); in the quasi-universal Gaussian 
model, this total energy is given by ~ e 9^ (ZM). A summary 
of our results are as follows: 

• Varying the characteristic jet width, 9„: We try a model in 
which 9 in equation 2 varies as a log-normal with mean 
< \g{9 ) > and standard deviation <J\ % (e„)- In general, we 
find that < lg(t9 /l°) > in the range 0.7-0.8 and cr\ %{ e„) in 
the range 0.3-0.5 does the best job reproducing the data. 
A higher mean tends to create an unobserved flattening 
at low values of 9j, while a lower mean causes the sim- 
ulated data to be too steep (with power law index much 
steeper than -2). Figure 4 shows the simulated data vs. 
actual data in the BFK range for < lg(t9 /l°) >= 0.8 and 
cig(#„) = 0.4. The simulated data shown in Figure 4 are fit 
with a power-law of index -1 .8 ± 0.2 in the range of the 
BFK data; this is consistent with the index we obtained 
when fitting a simple-power-law to the BFK data. The 
fit gives a reduced x 2 = 4.0, which - although formally 
unacceptable - is consistent with the value of x 2 obtained 
when we fit a power-law directly to the BFK data. Hence, 
we consider this similar value of k and x 2 to indicate that 
we are accurately reproducing the slope and scatter of 
the data. Finally, when we perform a KS test on the sim- 
ulated and observed data, we find that the probability that 
these data are derived from the same distribution is 69 %. 
A very low value of this probability << 1% would in- 
dicate with statistical robustness that these data are from 
different distributions. Hence, we consider this model to 
adequately reproduce the data. 

i and quantitatively similar results if we draw our 9 values from, say, a uniform 
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Using a normal distribution to produce variation in 6 , 
with average value < 9„ > and a standard deviation <jq , 
we find that values < 9 >« 7- 10° and a e „ ~ 2-4° do 
the best job in reproducing the data. Our best fit was 
for < 9 >= 9° and ag o = 4°. The simulated data in this 
case are fit with a power-law of index -2.0 ± 0.2 in the 
range of the BFK data; this is consistent with the index 
we obtained when fitting a simple-power-law to the BFK 
data. The fit gives a reduced x 2 = 6.4. A KS test gives 
a probability of 15 % that these data are derived from 
the same distribution. In general, this particular model is 
quite unstable (for small variations in 9 , we can get huge 
variations in the energy) and reproduces the data only in 
special cases. 

Varying the normalization, e : We vary the normalization 
e„ as according to a log-normal distribution with mean < 
lg(e G ) > and standard deviation eri g ( e „). For a 9„ = 10°, < 
lg(e o /10 50 erg) > in the range 3.4-3.7 and a cr lg(eo / 10 5o erg ) 
in the range 0.5-0.7 describes the data well. Figure 5 
shows the simulated data for < lg(e o /10 50 erg) >= 3.6 and 
^lgfe/io^erg) = 0-7 The simulated data shown in Figure 5 
are fit with a power-law of index -1 .5 ± 0.3 and x 2 = 5.6 
in the range of the BFK data. The formal power-law fit 
is marginally consistent with the fit to the BFK data, but 
a KS test performed on the simulated and observed data 
sets gives a probability of 99 % that the two data sets are 
derived from the same parent distribution. 2 

When we vary the normalization e according to a normal 
distribution with mean < e > and standard deviation a €o , 
we find that for a 9 Q = 10°, < e o /10 50 erg >= 1.3 x 10 3 
and cr €c = 1.7 x 10 3 describes the data well. Note that the 
values of the mean and dispersion cause some deviates to 
be negative. In this case, we take the absolute value of 
all deviates. The data here are fit with a power-law of in- 
dex -2.0 ±0.2 and x 2 = 7.4 in the range of the BFK data. 
A KS test performed on the simulated and observed data 
sets gives a probability of 49 % that the two data sets are 
derived from the same parent distribution. 

Varying 9„ and e„ in a correlated way: Finally, we vary 
9 and e„ such that the total energy e 9 1 = constant, 
again motivated by the result of Frail et al. (2001) and 
BFK that the total GRB energy is approximately con- 
stant. To do this, we take e from a log-normal distri- 
bution as described in the previous section and solve for 
9 from the above relationship. We find the observed 
data is reproduced when < lg(e o /10 50 erg) >= 3.4 and 
^lgtc/io^erg) = 0-6- This results in a mean jet width of 
20°. The simulated data shown in Figure 6 are fit with a 
power-law of index -1 .8 ± 0.2 and x 2 ^6.6 in the range 
of the BFK data. A KS test performed on the simulated 
and observed data sets gives a probability of 63 % that 
the two distributions are the same. 

find the observed 
>= 1.3 x 10 3 and 



For a normal variation in e a , we 
data is reproduced when < e o /10 50 



erg 

cr e„/io 50 erg = 2.0 x 10 3 . This gives a mean jet width of 
about 10°. These simulated data are fit with a power-law 
of index — 1.7 ± 0.2 and x 2 = 4.4 in the range of the BFK 
data. A KS test performed on the simulated and observed 



data sets gives a probability of 21 % that the two distri- 
butions are the same. 

The bottom line is that the data can be adequately reproduced 
with a Gaussian model for the jet structure, given some varia- 
tion in the models parameters (we note that in our simulations, 
we have assumed that the viewing angle is the agent that de- 
termines the jet break time). Each model above was able to 
reproduce the power-law behavior of the BFK data, the scatter 
(reflected by the value of x 2 of the power-law fit), and a KS 
test probability that shows that the simulated data are not in- 
consistent with being derived from the same parent distribution 
as the observed data. It appears, however, that the log-normal 
variations in e and the correlated e„ - 9 variations do the best 
job in reproducing the data (see Figures 5 and 6). The model 
in which 9 varies in general exhibits too much curvature and 
is less stable to the degree of variation. Although the correlated 
£o _ 9 variations do a good job of reprocing the data, the ad- 
ditional constraint that e„# 2 = constant may be considered an 
unattractive feature of this model. However, we do point out 
that this constraint is motivated by the observations (Frail et al., 
2001, BFK) and may in fact be a physical feature of GRBs. 

An interesting aspect of the Gaussian model is that the k = —2 
power-law behavior is not extended below values of 6j <J 2° and 
in some cases above 6j ;> 30°. Figure 7 shows simulations for 
the correlated e -9 variations, extendend above and below the 
range of the BFK data. A turnover is clearly seen at low values 
of 9j. For the case of either e„ or 9 varying individually, or 
for the case of the correlated e - 9 with a low value for the 
average jet width <^5°, we also find a cutoff at high values of 
9j :> 30°. We note that Figure 7 uses the same parameters as 
in Figure 6, where the average jet width is 20°. This average 
jet width can vary and this will affect exactly where the low an- 
gle turnover and high angle cutoff appear. Nonetheless, if this 
model were correct, then we should not see many more bursts 
with 9j <^ 1°, and possibly also 9j J> 30°. These predictions may 
be tested with future observations of GRBs, and in fact the low 
9j prediction may have already been seen. Interestingly, BFK 
discuss several outliers in their sample that appear to be "fast- 
fading" (possibly because the jet break occured at very early 
times, implying a very small value for the jet opening angle 6j) 
and sub-luminous relative to the rest of the sample. This is con- 
sistent with the trend predicted by the Gaussian jet structure at 
small values of 9j. 

4. LUMINOSITY FUNCTIONS 

4.1. Predictions 

One of the advantages to the quasi-universal jet configura- 
tion, is that it makes definite predictions for the observed lu- 
minosity function (LF) of GRBs. For example, as discussed in 
ZM, the isotropic luminosity function N(e) at a given redshift 
can be determined from the burst angular distribution N(9)d9 oc 
sin9d9 oc 9d9. In the power-law model, 9 oc e~ l l k , so that 
N(e) oc 9d9/de oc et" 1 " 2 /*). For a Gaussian model, N(e) oc e~ l . 
These predictions are derived in the case of no variations in the 
model parameters. We extend this analysis by simulating lumi- 
nosity functions, allowing parameters to vary in our power-law 
and Gaussian models. Using the best fit parameters and varia- 
tion schemes from the previous analysis, we simulate 500,000 



2 Again, we note that the KS diagnostic is intended to check whether two distributions can be considered statistically different. A value << 1 % would indicate 
such a result. A higher value (like 99 %) does not necessarily indicate the model is any better than a lower value (like 50 %). All of our KS probabilities show that it 
is not the case that the simulated and observed data are derived from different distributions. 
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GRBs for each jet structure with the distribution of observing 
angle following N(6 v )d6 v cx sm9 v dd v , and where the range of 
the observing angle is from 0° to 90°. Again, these simulations 
hold for luminosity functions at a given redshift. 

4.2. Power-law Jet Structure 

We first test our methods by simulating the GRB LF in the 
case where e(8) is a simple power-law (with no parameter varia- 
tions, and our best fit index of k = -1 .9). As expected we obtain 
a slope of the LF of index w -2 (solid line in Figure [5}, con- 
sistent with the analytical expression of ZM and given above 
in §4.1. The minimum luminosity of ~ 10 51 erg is due to the 
upper boundary of the observer angle 8 V = 90°. We then simu- 
late the cases where the power-law index k varies according to 
a normal distribution, and e„ varies according to both normal 
and log-normal distributions. The results are shown as dotted, 
dashed, and dash-dotted lines in Figure[8] respectively. The lu- 
minosity functions obtained when we allow parameters to vary 
are similar to the simple power-law jet at the high luminos- 
ity end, being well described by a power-law of index -2 for 
more than three order of magnitudes. At the low luminosity 
end, however, the LFs tend to flatten. This flattening is a result 
of the low-luminosity cutoff present in the simple power-law 
case. For example, above the cutoff, when we allow parame- 
ters to vary there is a compensation between lower and higher 
luminosities that preserves the overall slope of the luminosity 
function. However, variations in the parameters near the cutoff 
cause some GRB luminosities to be pushed below the mini- 
mum luminosity, without compensation from GRBs of lower 
luminosity. 

4.3. Gaussian Jet Structure 

Again, we first simulate the luminosity function of a Gaus- 
sian jet, with fixed normalization, e , and characteristic jet 
width, o . The result is a powerlaw of slope -1, shown by the 
solid line in Figure|9] and consistent with the analystical expres- 
sion in §4. 1 . Note that this luminosity function has a maximum 
luminosity corresponding to an observer angle of 0°. We then 
vary the characteristic angle, 6 , of the Gaussian jet according 
to a log-normal distribution. The simulated luminosity func- 
tion is very similar to the luminosity function obtained from 
constant Gaussian jet model; it has a power-law index of -1 
and an upper bound, shown as dotted line in Figure|9] Next, we 
consider the case where e varies as log-normal. The result is 
shown by the dashed line in Figure|5] The luminosity function 
resembles the case of constant e Q and o on the low luminos- 
ity end, but steepens at high luminosities. This steepening is 
due to the same type of effect that caused the flattening at the 
low luminosity end in the power-law jet case (see §4.2 above) 
- variations in the parameters near the maximum luminosity 
cause some GRB luminosities to be pushed above the maxi- 
mum, without compensation from GRBs of higher luminosity. 
We also simulate the luminosity function in the case where both 
e and o vary, but preserve the relation e a 8l = constant (see §3). 
The luminosity function is shown as dash-dotted line in Fig- 
ure|5] The LF in this case has a slightly flatter power-law index 
at low luminosities, with a clear break at L ~ 10 52 erg and a 
sharp steepening after that. This break occurs because an anti- 
correlated e„ - 6 distribution causes an enhanced pile-up near 
(< lg(e G ) >, < lg(0 o ) >) in the e— 6 plane, leading to a flattening 
of the LF at the corresponding luminosity. Finally, we simulate 
LFs from Gaussian jets where the parameters vary according 
to a normal distribution; the LFs are shown in Figure Q75| The 



qualitative behavior is similar to the log-normal variation case, 
although these LFs are steeper at the high luminosity end over 
a smaller dynamic range. 

The luminosity functions in both the power-law and Gaus- 
sian jet structures have a definite change of power-law index 
when variations in the model parameters are included. Power- 
law jets produce LFs with a slope of -2, but with a dramatic 
flattening at luminosities below about 10 51 erg. Gaussian jets 
produce LFs with a slope of <~ — 1, although varying e , and 8 
and e a in an anti-correlated way, causes the LF to steepen at 
high luminosities. In particular, the anti-correlated e -8 case 
with log-normal parameter variations produces a clear break at 
<J 10 52 erg, above which there is a steepening with an index of 
~-2. 

4.4. Comparing with observed data 

We would like to compare the predictions for the GRB lu- 
minosity function from the quasi-universal jet model of GRBs, 
with the observed GRB luminosity function. Unfortunately, we 
are not yet at the point where we can directly measure the GRB 
LF with any statistical certainty. This is because there are still 
on a handful (~ 33) of GRBs with measured redshifts, and de- 
termination of the LF requires large numbers (at least 100) of 
luminosities from a uniform sample (and a sample with com- 
plete and detailed knowledge of the selection effects). 

Usually, the GRB LF is estimated by utilizing a large sam- 
ple of GRBs with a measured flux distribution n(f), assuming 
an underlying density distribution, p(z), cosmological model, 
and then computing the luminosity function N{L) through the 
equation n(f)df = N(L)p(z)/(l+z)(dV /dz)dLdz. Note that this 
assumes independence of the variables L and z, or - in other 
words - that the luminosity function does not evolve with red- 
shift. Several studies have taken this approach to estimate the 
GRB LF. For example, Stern et al. (2003) - assuming a steep 
decline in the GRB source population above a redshift of 1.5 - 
find N(L) oc ZT 1 - 4 up to L ~ 3 x 10 51 erg and then sharply de- 
clines after that. Schmidt (2002) - using various models for 
the GRB rate density based on the estimates for the global star 
formation rate - find that the GRB luminosity function can be 
described by a broken power-law with an index ~-l on the 
low luminosity end (< 3 x 10 51 erg) and sharply declines after 
that, consistent with the Stern et al. results. [All luminosi- 
ties quoted are the isotropic equivalent luminosities.] Although 
these results are fairly uncertain because they rely on the as- 
sumption of the underlying GRB rate density, they appear to 
be marginally consistent with predictions from the Gaussian jet 
structure. The power-law universal jet may also be consistent 
with these results if the data is fit with a single power-law such 
that the steepening of the LF above some break causes an over- 
all steepening of the LF index, closer to -2. [We note that re- 
cently, some attempts have been made to constrain the GRB LF 
based on the small sample of bursts with direct spectroscopic 
redshifts. For example, van Putten & Regimbau (2003) find 
that a log-normal description of the LF does a suitable job of 
describing the existing data.] 

Although spectroscopic redshifts have only been obtained to 
a number of GRBs too small to determine a luminosity function 
directly, in recent years there has been evidence of so-called 
"luminosity indicators" (Fenimore & Ramirez-Ruiz, 2000, Re- 
ichart et al., 2001, Norris et al. 2001) from which luminosities 
of GRBs can be obtained. For example, Fenimore & Ramirez- 
Ruiz (2000) found a correlation between GRB luminosity and 



6 



the light curve variability (an estimate of the "spikiness" of the 
light curve) based on a small sample of bursts with measured 
redshifts. From this, they were able to obtain luminosities and 
redshifts to 220 GRBs. Because they used BATSE gamma-ray 
data and chose a flux limit to their sample, the selection ef- 
fects were very well defined and therefore a luminosity function 
can be estimated. Fenimore & Ramirez-Ruiz assumed indepen- 
dence between L and z in their sample and found a luminosity 
function with a power-law index of approximately -2, consis- 
tent with the power-law quasi-universal jet configuration. A 
more detailed analysis of this sample, however, shows that L 
and z are in fact correlated (Lloyd-Ronning, Fryer, & Ramirez- 
Ruiz, 2002, Yonetoku et al. 2003, Graziani et al. 2003); this 
correlation has to be accounted for to correctly determine a 
LF from the data. Lloyd-Ronning et al. did this by defining 
a variable U = L/X(z), where A(z) parameterizes the correla- 
tion between L and z- They found a single power-law fit gives 
N(L') oc L'" 2 - 2 , while a broken power-law fit gives N(L') oc L'~ l 5 
for U < L' a and N(L') oc L'" 3 3 for L' > L' a , where L' = 5 x 10 51 
erg. Again, there is no underlying assumption about the density 
distribution of GRBs in this analysis and the correlation be- 
tween L and z present in the data has been accounted for. [We 
mention that in the context of the quasi-universal jet model, 
luminosity evolution could imply the overall normalization e c 
evolves with redshift, or that the characteristic opening of the 
jet o evolves.] The single power-law result appears to be 
more consistent with the power-law quasi-universal jet struc- 
ture, while the broken power-law result is better described by 
the (at least qualitative) predictions of the Gaussian case. Fi- 
nally, Norris (2002) found - using the luminosities and redshifts 
obtained from the luminosity-lag relation, another luminosity 
indicator (Norris et al, 2001) - that the GRB LF scales as L~ l 
for low luminosities and L~ l s for high luminosities. This is 
very similar to the results we obtain for the LF from a Gaussian 
jet when e a and o vary in an anti-correlated way. 

Although some of the above results suggest the plausibility 
of a Gaussian jet structure for GRBs, the predictions of the 
quasi-universal jet model GRB LFs are better tested when more 
GRBs with redshifts are obtained and a direct GRB LF can be 
measured. This is a primary goal of the Swift satellite, to be 
launched in June of 2004. 

5. SUMMARY AND DISCUSSION 

In this paper, we have tested two possible models for a quasi- 
universal jet structure to a GRB. In particular, we attempt to 
reproduce the data of BFK who showed an anti-correlation be- 
tween the isotropic emitted energy Ei so and the jet opening an- 
gle, Oj derived in a uniform jet model, such that Ej SO oc 9~ 2 . 
In the quasi-universal jet model (Rossi et al., 2002, ZM), in 
which the jet opening angle and total emitted energy are ap- 
proximately the same for all GRBs, but the energy varies as 
a function of angle from the jet axis within a GRB, this cor- 
relation is a reflection of the energy profile e(8) as a function 
of angle from the jet axis. Motivated by ZM, who suggested 
that this profile may take on different functional forms, we have 
tested how Gaussian and power-law models for e{6) reproduce 
the BFK data. In particular, we have allowed for realistic scat- 
ter in the parameters of each model (letting each parameter vary 
as a either a log-normal or normal distribution with some mean 
and standard deviation). 

We find that both the power-law and Gaussian jet structures 
can reproduce the data quite well, with only minimal scatter 



required in the model parameters - particularly when the pa- 
rameters vary according to the more physically reasonable log- 
normal distribution. The strengths of the power-law model are 
its simplicity and ease in reproducing the observed data. How- 
ever, as discussed in the introduction and throughout the text, 
there are several reasons to consider other jet configurations, 
and in particular the Gaussian model appears to adequately de- 
scribe the observations. Furthermore, The Gaussian model pre- 
dicts a sharp turnover in the E iso - 9j plane for low (and in some 
cases high) values of Oj. This turnover at low jet angles has pos- 
sibly been observed by BFK, who found a few "sub-luminous" 
bursts (relative to the rest of their sample), with steeply de- 
clining light curves which could indicate a very small opening 
angle 9j. These "outliers" are consistent with the trend pre- 
dicted by the Gaussian quasi-universal jet configuration at low 
8j. Furthermore, these outlier bursts challenge the BFK con- 
clusion that the energy resevoir in GRBs is approximately con- 
stant, when taken in the context of the uniform jet model. How- 
ever, the measured Oj is - in the context of the quasi-universal 
jet paradigm - the viewing angle and not the characteristic width 
of the jet 9 , which could be larger. This means that the total 
energy within the GRB could still be standard, if the results are 
considered in the quasi-universal jet picture with a Gaussian jet 
structure. We also find that the luminosity function predicted by 
the Gaussian model, including variations in the model parame- 
ters, appears to be consistent with past studies of the luminosity 
GRB LF, although these predictions are better tested when the 
GRB LF can be directly measured. 

Perna et al. (2003) showed that a universal jet model with a 
power-law structure (of index of -2) predicts a distribution of 
break times in the afterglow light curve that is consistent with 
what is observed. It would be interesting to explore their analy- 
sis in the context of the Gaussian model, which is proving to be 
a viable model for the quasi-unviersal jet structure of GRBs. 
(Given our results above, we suspect that a Gaussian model 
may give qualitatively similar results as the power-law.) We 
note that Lamb et al. (2003) have recently pointed out that the 
quasi-universal jet model fails to reproduce the large dynamic 
range of the observed relationship between isotropic energy and 
spectral peak energy (see, e.g., Amati et al., 2002), which spans 
not only the "classic" GRB energies (i.e., those of BATSE, from 
~ 50 keV to 1 MeV), but also includes so-called X-ray flashes, 
which have spectral peak energies down to a few keV. Their 
conclusions are in the framework of a power-law structure for 
the quasi-universal jet, and such an investigation in the frame- 
work of the Gaussian model is underway (Zhang et al., in prep). 
Furthermore, we comment that the realistic jet struture might 
not be strictly power-law, Gaussian, nor their simple superpo- 
sitions; in fact, recent numerical simulations (Zhang, Woosley, 
& MacFadyen, 2003) indicate a double Gaussian structure for 
the jet may provide the best description (where one Gaussian 
is used for the core of the jet and one for the wings). The bot- 
tom line is now that we can reproduce the data with varying 
parameters in these simplified models, we can do it with more 
realistic models too. And such models may be able to explain 
all of the observed gamma-ray burst data in the context of a 
quasi-universal jet configuration. 

We would like to thank Dale Frail for sending us a table of 
the BFK data. We would also like to thank the referee for in- 
sightful comments and suggestions that improved this work. 
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FIG. 2. — Isotropic emitted energy vs. jet opening angle. The square points are the BFKdata, while the circles indicate data derived from a power-law jet structure, 
with an index k that varies according to a normal distribution with mean of -2 and standard deviation of 0.7. 
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FIG. 3. — Same as Figure 2, but where the circles indicate data derived from a power-law jet structure, where the normalization varies as a log-normal distribution, 
with a mean < lg(e„/10 50 erg) >= 4.8 and a standard deviation &\ g ( to H(pa ag) = () f >- 



10 4 



- o 



_ o 



fcooo 



100 



o 



o, 



c 



o 



o 



o 



■ -i 



o 



o 



c 



10 

10 

B } (degrees) 

FIG. 4. — Same as Figure 2, but where the circles are data derived from a Gaussian jet structure, with a characteristic width O that varies as a log-normal 
distribution with a mean < lg(6 /l ) >= 0.8 and a standard deviation fjg^/io) = 0.4. 
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FIG. 5. — Same as Figure 2, but where the circles are data derived from a Gaussian jet structure, with a characteristic width 9„ = 10° and a normalization e„ that 
varies as a log-normal with a mean < lg(e o /10 50 erg) >= 3.6 and a standard deviation Ojg^/iosOag) = 
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FIG. 6. — Same as Figure 2, but where the circles are data derived from a Gaussian jet structure, a normalization e that varies as a log-normal with a mean 
< lg(e o /10 50 erg) >= 3.4 and a standard deviation Ci g ( eo /io5o erg ) = 0.6, and a characteristic jet width that follows the relation e o 0g = constant. 
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FIG. 7. — Same as Figure 6, but plotted over a larger range in both axes to show the turnover at low 9j. For the case of < e > or < O > varying individually, 
there is also a sharp cutoff at high 6j > 30 degrees. 
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FIG. 8. — Monte-Carlo simulated (differential) luminosity functions for GRBs with power-law jets, where No is an arbitrary normalization. The luminosity 
function is calculated assuming isotropic emission and a typical duration to all bursts. The solid, dotted, dashed, and dash-dotted lines represent power-law jets with 
no variations, with powerlaw index k varying as according to a normal distribution, with eo varying as according to a normal distribution, and with eo varying as a 
log-normal, respectively. 
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FIG. 9. — Same as Figure 8, but for Gaussian jets and log-normal parameter variations. The solid, dotted, dashed, and dash-dotted lines represent Gaussian jets 
with no variations, with 9q varying, with eg varying, eg and 9q varying in an anti-correlated way, respectively. 
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FIG. 10. — Same as Figured but with the parameters varying according to a normal distribution. 



